RevMexAA (Serie de Conferencias), 00, [l]-g (2000) 



SPECTROPHOTOMETRY OF MASSIVE STAR FORMING REGIONS WITH 

GTC 

M. Castellanos and A. I. Diaz 
Universidad Autonoma de Madrid, Spain 
and 
E. Terlevich 1 
INAOE, Mexico 

RESUMEN 

En estos ultimos afios, la detection de estrellas Wolf-Rayet (WR) en regiones HII gigantes ha levantado ciertos 
interrogantes sobre nuestros conocimientos de la evolution de estrellas masivas y sus atmosfcras extensas, la 
edad de la poblacion ionizante y su impacto sobre las propiedades ffsicas de dichas regiones. En este arti'culo, 
presentamos observaciones espectrofotometricas de 4 regiones HII extragalacticas gigantes que muestran rasgos 
debidos a estrellas WR en sus espectros de emision. Nuestro objetivo consiste en reproducir de forma simultanea 
tanto las propiedades de las estrellas WR observadas, como el espectro de lfneas de emision del gas, con la 
ayuda de modelos de smtesis de poblaciones. Finalmente, senalamos las principales ventajas que proporcionara 
GTC para nuestro mejor entendimiento de regiones de formation de estrellas masivas. 

ABSTRACT 

In these latest years, the detection of Wolf-Rayet (WR) stars in Giant HII Regions (GHRs) has yielded several 
questions about our current understanding of massive stars evolution and hot expanding atmospheres, the 
age of the ionizing populations and their impact onto the physical properties of GHRs. Here, we present 
spectrophotometric observations of 4 extragalactic GHRs which show WR features in their spectra. Our goal is 
to reproduce simultaneously the observed WR properties and the emission line spectra with the help of current 
evolutionary synthesis models. Finally, we address the main advantages that GTC will provide to our better 
understanding of massive star forming regions. 
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1. INTRODUCTION 

WR stars are evolved descendants of massive O 
stars, and hence, they represent one of the latest 
stages of massive star evolution. These stars exper- 
iment powerful winds, which can lead to a complete 
loss of their outer envelopes. Let us summarize the 
main topics associated with the presence of WR stars 
in GHRs. Stellar evolution predicts mass loss rates 
to be more prominent at high metallicity (Z > Z ) 
and, therefore, O stars can enter the WR phase at 
a lower cutoff mass (Meynet 1995). This leads to 
both higher WR/O star ratios and stronger recom- 
bination lines formed in the wind. However, few de- 
tections at high metallicity have been made at the 
moment (Castellanos, Diaz & Terlevich 2002, Paper 
I; Bresolin & Kennicutt 2002) to shed some light on 
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this matter. On the other hand, a high metal content 
implies an increase in the opacity of the stellar ma- 
terial. This effect would lower the effective tempera- 
ture of massive stars in regions of high metal content. 
There is a general agreement about the hardening of 
the ionising radiation in regions of low metal con- 
tent (Campbell et al. 1986). However, the soften- 
ing of the ionising radiation in regions of high metal 
content is difficult to quantify, since the functional 
parameters in these regions cannot be parametrised 
in a trivial way (see Diaz et al. 1991). Stellar evo- 
lution can also predict the WR/O star ratio from 
the elemental surface abundances and, therefore, it 
should be able to reproduce adequatelly the WR 
population properties (mass loss rates, strength of 
the lines). Again, the detailed observation of GHRs 
with embedded WR stars is compulsory in order to 
constrain the stellar evolution assumptions. Finally, 



2 



CASTELLANOS ET AL. 



taking apart the metal content of the regions, the 
sole presence of WR stars is supposed to rise up 
drastically the number of ionising photons at ener- 
gies higher than ~ 40 eV (Perez 1997). Hence, one 
might expect the ionization structure to change in- 
side GHRs. In order to analyze these topics, we have 
studied 4 GHRs in NGC 628 (H13) and NGC 1232 
(CDT1, CDT3 and CDT4) showing WR features in 
their spectra (Castellanos, Diaz & Terlevich 2002, 
Paper II). 

2. WR POPULATIONS 

Schaerer & Vacca (1998; SV98) have presented 
very detailed models of the WR population in young 
star clusters, at different metallicities from Z=0.001 
(1/20 Z Q ) to Z=0.04 (2 Z ). Their clusters are 
formed according to a Salperter IMF with upper and 
lower mass limits of 120 M Q and 0.8 M Q respectively. 
They use the stellar evolution models by Meynet ct 
al. (1994) which assume a mild overshooting and 
enhanced mass loss rate. These models have been 
shown to reproduce the observed WR/ O star ratios 
in a variety of regions (Maeder & Meynet 1994). Re- 
garding the energy output of main sequence stars 
they use the CoStar models (Schaerer & de Koter 
1997) which include non-LTE effects, line blanket- 
ing and stellar winds, for stars with initial masses 
larger than 20 M Q and Kurucz (1992) plane-paralel 
LTE models, including line blanketing effects, for less 
massive stars. The atmospheres of evolved stars in 
the WR phase correspond to the spherically expand- 
ing, non-LTE, unblankcted models by Schmutz, Lei- 
therer & Gruenwald (1992). 

SV98 provide accurate predictions, as a func- 
tion of the cluster age, for the total number of WR 
stars and their subtype distribution, the broad stel- 
lar emission lines and the luminosities and equiva- 
lent widths of the two "WR bumps" at A 4650 and 
A 5808 A . We have used these models to derive the 
age of the ionising population which contains WR 
stars and whose metallicity has been previously de- 
rived by means of the calculated ionic temperatures 
(see Paper I). 

Figure 1 shows the predicted emission line inten- 
sities and equivalent widths of the WR blue "bump" 
for different metallicities (0.2Z Q ; 0.4Z© and Z Q ) to- 
gether with the observed values in the four regions. 
Solid symbols correspond to the observed data. It 
can be seen that in three of the analyzed regions, 
observations and predictions are in excellent agree- 
ment for a single instantaneous burst. The derived 
ages for the WR population range between 3.1 Myr 
(CDT3) and 4.1 Myr (H13 and CDT4). In fact, for 



regions CDT3 and CDT4 observations and predic- 
tions arc nearly identical for the respective clusters 
and therefore only one symbol is shown. In the case 
of region CDT1, two open symbols are shown, one 
corresponding to the predictions for a single ionising 
cluster of 2.4 Myr, and another one including the 
contribution to the H/3 and continuum luminosity of 
a cluster not containing WR stars (around 7 Myr). 
This latter one is able to reproduce the observed WR 
features and the H/3 equivalent width. 

3. MODELLING THE EMISSION LINE 
SPECTRA 

Regarding the integrated Spectral Energy Dis- 
tribution (SED) from the ionising cluster, the recent 
models from Leitherer et al. (1999; STARBURST99) 
provide an almost sclfconsistcnt frame to be used in 
combination with the WR models described above. 
In fact, they use the same stellar evolution mod- 
els with enhanced mass loss rate, the same atmo- 
sphere models to describe the stars in the WR phase 
and cover the range of metallicities and IMF used 
by SV98. The only appreciable difference between 
both sets of models concerns the atmospheres of the 
main sequence stars with initial masses greater than 
20 M© which in STARBURST99 are represented by 
the plane-paralel Kurucz 's models implemented by 
Lejeune, Cuisiner & Buser (1997). We have there- 
fore used the STARBUST99 models in order to fit 
the emission line spectra for our analysed regions. 
Hence, as the WR emission features have been used 
to date the ionizing star clusters, these should repro- 
duce the observed gas emission line intensities if the 
other parameters controlling the emission line spec- 
tra, namely elemental abundances, particle density 
and ionisation parameter, are known. Our four ob- 
served GHRs meet all these requirements. 

Therefore, since both SV98 and STARBURST99 
models make use of the same stellar evolution pre- 
scriptions, we have assumed that a STARBURST99 
model of a cluster of a given metallicity which con- 
tains the same WR/O star ratio as SV98 (corre- 
sponding to a given age), must provide the spectral 
energy distribution of the ionising radiation. In the 
case of region H13 in NGC 628, this model (corre- 
sponding to an age of 4.1 Myr from the previous 
analysis of the WR features) however provides an 
[OIII]A5007 A line emission which is higher than ob- 
served by a factor of about 4. On the other hand, our 
best photoionisation model corresponds to an age of 
4.7 Myr producing an H/3 equivalent width of 108 A, 
to be compared with the observed one of 140 A. The 
rest of the spectrum is reproduced remarkably well 
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Fig. 2. Relative intensity (left) and equivalent width (right) of the WR blue 'bump' versus H/3 equivalent width for 
SV98 models of three different metallicities as labelled. The data are shown as solid and open symbols as explained in 
the text. 



taking into account the uncertainties in both the ob- 
served values and model computations. We can con- 
clude that a single star cluster with age between 4.0 
and 4.7 Myr fit all the observations adequately. 

For region CDT1, given that a young 2.4 Myr 
cluster is successful in reproducing both the WR blue 
'bump' luminosity and equivalent width (although 
overpredicts the equivalent width of H/3 by a fac- 
tor of seven), a composite population with, at least 
two clusters must be invoked. We have therefore run 
a model using as ionising source the combination of 
the spectral energy distributions of two ionising clus- 
ters of 2.4 and 7.1 Myr of age calculated with the 
STARBURST99 code in which the youngest of the 
two provides 10 times the number of ionising photons 
emitted by the oldest. This cluster contains the same 
WR/O star ratio as the SV98 model reproducing the 
observed WR features. The computed photoionisa- 
tion model is able to reproduce all the observables 
adequately. 

For the other two GHRs in NGC 1232, CDT3 
and CDT4, STARBURST99 model clusters with the 
same WR/O ratio as those given by SV98, provide 
[OIII]A5007 A line emission higher than observed 
by factors of about 5 and 3 respectively. Clusters 
slightly older, also reproducing the observed WR 
features within the errors, have larger WR/O num- 
ber star ratios (up to 8.3 xl0~ 2 ) and therefore their 
spectral energy distributions result even harder. Al- 
ternatively, it is possible to find a combination of 
ionising clusters which reproduces well the emission 
line spectrum, but predicts WR feature luminosities 
and equivalent widths well below the observed val- 
ues. 



4. DISCUSSION AND CONCLUSIONS 

The fact that the WR features are adequately 
reproduced by SV98 models seem to imply that the 
evolutionary tracks are able to predict the right rel- 
ative numbers of WR/O stars, and their different 
subtypes at the derived abundances. These relative 
numbers, combined with the observed emission line 
luminosities of the individual WR stars, and the pre- 
dicted continuum energy distribution of the ionis- 
ing population predict emission line intensities and 
equivalent widths of the WR stars that are in ex- 
cellent agreement with observations. On the other 
hand, in two of the analyzed regions (CDT3 and 
CDT4 in NGC 1232), it is not possible to fit the 
emission line spectrum since the population contain- 
ing WR stars produces a spectral energy distribution 
which is too hard to explain the emission of the gas. 
The same sort of effect has been found by Diaz et 
al. (2000a) for region 74C in NGC 4258 and Este- 
ban et al. (1993) for the galactic WR nebula Ml-67. 
These latter authors, from stellar and nebular spec- 
troscopic analyses, concluded that lower tempera- 
tures were required from the photoionisation models 
for late type WN (WNL) stars. In their study they 
used the unblankcted WR models of Schmutz et al. 
(1992). A subsequent reanalysis of this region has 
been made by Crowther et al. (1999) using blan- 
keted model atmospheres. In this case the result- 
ing ionising spectrum is much softer and reproduces 
better the observations although some discrepancies 
still remain. For the high metallicity GHR CDT1 
in NGC 1232, a composite population can explain 
adequately both the WR features and the emission 
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line spectrum. Composite populations for HII re- 
gions have been found in previous works by Mayya 
& Prabhu (1996) and by Diaz et al. (2000b) for disc 
and circumnuclear objects respectively from broad 
band and Ha photometry. 

In the case of region H13 in NGC 628, a single in- 
stantaneous burst between 4.0 and 4.7 Myr, is able to 
reproduce all the observables within the errors. This 
result seems to indicate that line blanketing effects at 
low metallicity (0.2 Z©) could be less severe for the 
correct interpretation of the emission line spectra. 

Our observations indicate that no appreciable 
change is found in the ionization structure of the 
analyzed HII regions, despite the presence of WR 
stars. Line blanketing in WR atmospheres would 
point again in the right direction. It should be kept 
in mind that all our analysis is based on ionisation 
bounded models. An approach of the emission line 
spectra of HII region with WR features on the ba- 
sis of matter bounded models has been presented in 
Castellanos, Diaz & Tenorio-Tagle (2002). Models 
of this kind have been found to provide excellent fit- 
tings to the observations of H13 in NGC 628, CDT3 
in NGC 1232 and 74C in NGC 4258 using both SV98 
and STARBURST99 models. These models would 
point to an important leakage of ionising photons 
depending on both the metallicity and evolutionary 
state of the region. There are still several caveats 
about our understanding of emission line spectra of 
GHRs with WR stars. Certainly, the controversial 
observation of main sequence hydrogen-rich WNL- 
like stars (de Koter, Heap & Hubeny 1997; Crowther 
& Dessart 1998), not predicted by current stellar evo- 
lution calculations at ages less than 2Myr at low 
metallicity, and the observation of dust envelopes 
around WC stars (Marchenko et al. 2002) may have 
also important consequences to interpret correctly 
the ionization structure in GHRs, closely related to 
the impact of WR stars in the surrounding medium. 
The combination of GTC + OSIRIS spectroscopy 
would allow the observation of WR features in low 
excitation HII regions (high metallicity regions, in 
general) at shorter exposure times. The detection of 
He II at 4686A with Tunable Filters seems to be a 
promising tool with GTC. Clearly, GTC would allow 
to solve some of the problems mentioned above. 
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